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ABSTRACT 


We study the effects of neutron captures in AGB stars on “Fe-group” ele¬ 
ments, with an emphasis on Cr, Fe, and Ni. These elements show anomalies in 
54 Cr, 58 Fe, and 64 Ni in solar-system materials, which are commonly attributed 
to SNe. However, as large fractions of the interstellar medium (ISM) were re¬ 
processed in AGB stars, these elements were reprocessed, too. We calculate the 
effects of such reprocessing on Cr, Fe, and Ni through 1.5M 0 and 3M 0 AGB 
models, adopting solar and 1/3 solar metallicities. All cases produce excesses of 
54 Cr, 58 Fe, and 64 Ni, while the other isotopes are little altered; hence, the obser¬ 
vations may be explained by AGB processing. The results are robust and not 
dependent on the detailed initial isotopic composition. Consequences for other 
“Fe group” elements are then explored. They include 50 Ti excesses, and some 
production of 46 ’ 47,49 Ti. In many circumstellar condensates, Ti quantitatively 
reflects these effects of AGB neutron captures. Scatter in the data results from 
small variations (granularity) in the isotopic composition of the local ISM. For 
Si, the main effects are instead due to variations in the local ISM from differ¬ 
ent SNe sources. The problem of Ca is discussed, particularly with regard to 
48 Ca. The measured data are usually represented assuming terrestrial values for 
42 Ca/ 44 Ca. Materials processed in AGB stars or sources with variable initial 
42 Ca/ 44 Ca ratios can give apparent 48 Ca excesses/deficiencies, attributed to SNe. 
The broader issue of Galactic Chemical Evolution is also discussed in view of the 
isotopic granularity in the ISM. 

Subject headings: stars: AGB and post-AGB - Solar system: meteorites: Solar system: 
presolar grains - ISM: abundances - Nuclear reactions, nucleosynthesis, abundances 
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Introduction 


The purpose of this study is to investigate the effects of slow neutron captures in 
AGB stars on isotopes of the “Fe group” elements. In particular Ti, Cr, Fe, Ni, and Zn 


are considered to be members of the “iron group”. These nuclei are predominantly t 


product of supernovae 

(SNe), 

x)th core collapse and SNela 

see e 

Chieffi & Limongi 

2013; 

Langer 

2013; 

Hillebrandt & Niemeyer 

2000 


le 


Moll ct al. 


2014; 


and references therein). 


It has been known for over three decades that small (< 1%) variations are found in the 
isot opic compositi o ns of Ti, Cr, Fe, Ni, and Zn in macroscopic samp l es of m eteorites: see 


e.g. 


Ireland ([1990); 


Ireland & Fegley (120001 ) and the review by 


Birck (119 86). In particular, 


these isotopic effects are common in calcium-aluminum-rich inclusions (hereafter CAIs) in 
meteorites, but are also wide spread at lower levels in “bulk” samples of different groups 
of meteorites. The observed effects in Cr, Fe, and Ni typically show enrichments in the 
heaviest isotopes 54 Cr, 58 Fe, and 64 Ni, which are a small fraction of the abundance of the 
respective element (~ 2.3% for 54 Cr, ~ 0.28% for 58 Fe, ~ 0.9% for 64 Ni). In addition, the 
isotopic patterns observed in Ti, Ca, and Si are distinct and complex. Only one datum 
exists for Zn. As the Fe group elements are associated with Supernovae (SNe), it has been 
natural that interpretations of these effects in CAIs and “bulk” meteorites have focused 
on nucleosynthetic processes in them. This has not often led to any firm explanation of 
the observed isotopic anomalies or to clear predictions. Th e first s tudy that considered the 


possibility of s-process effects on Fe-group nuclei was by 


The et al. 


02Oo 3) , who explored the 


occurrence of s-processing in massive stars after revisions in the rate for the 22 Ne(cpn) 25 Mg 
reaction and provided some possible insights. 


The s-processed material is, in the interstellar medium (hereafter ISM), predominantly 
from low mass AGB stars. We expect that the effect of s-processing on Cr, Fe, and Ni 
(originally synthesized in SNe) is to produce the heavier isotopes of these elements, 54 Cr, 
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58 Fe, and 64 Ni, in over-abundance; this can be seen by considering that the pure s-process 
abundance pattern is, in general, smoothly correlated with the neutron-capture cross 
sections, yielding processed material enriched in nuclides heavier than the seed nuclei. 

As only a very small fraction of matter in the star is subjected to s-processing, it is not 
obvious just how large the consequences will be in the envelope of an AGB star, or will be 
provided to the ISM. There are thus three distinct issues: i) What are the s-process effects 
on isotopic abundance patterns in Cr, Fe, and Ni in an AGB envelope? ii) How significant 
are these effects? iii) What are the effects to be expected on the nuclei of lower atomic 
number than iron (e.g. Ti, Si, and Ca)? 

We will show that there are substantial consequences on the abundances of 54 Cr, 58 Fe, 
and 64 Ni produced in AGB stars and that these may provide an explanation of the effects 
observed in many macroscopic meteorite samples. We will try to relate this to data in 
presolar dust grains. The shifts found in other Fe-group nuclei subject to the s-process will 
be shown to be more complex. 

In the different classes of data that are available we must consider both presolar grains, 
each of which come from a particular stellar source, and macroscopic samples of solar 
system material, that come from blends from different stellar sources. 


The macroscopic samples of very early solar system matter derive from mixing, heating, 
melting from various precursor materials and exhibit relatively small isotopic shifts (< 1%). 
They are blends of very local parcels from unknown sources, but mixed so that they are 
rather close to the average isotopic admixture. In contrast to such macroscopic samples, 
grains formed as circumstellar condensates (subsequently called CIRCONs) and found in 
meteorites may exhibit very large (> 50% or more) isotopic abundance shifts (see e.g. 


Z i nner 


20141 for an in-depth review). We will try to relate the results on Cr, Fe, and Ni to 


data on presolar dust grains. A CIRCON is the direct product of stellar nucleosynthesis in 





5 


a single star. The initial material from the ISM that is incorporated in a forming star may 
be variable, but mixing provides an initially homogeneous isotopic pattern in that star, 
which is subject to processing in the subsequent stellar evolution before the formation of 
a grain. The effects produced will depend on the particular stage in the evolution and on 
the initial isotopic composition of the star. Results on CIRCONs often serve as clues to 
still unrecognized stellar processes. There are often difficulties in explaining nucleosynthesis 
effects in a circumstcllar dust particle from a single star resulting from incomplete knowledge 
of many nuclear and dynamical processes. For an individual CIRCON grain, the isotopic 
composition then depends on the initial abundances (not known) of the associated star and 
on the details of the subsequent stellar evolution which we claim to understand. For stellar 
sources that have, at each time, a well mixed envelope, the problem is simplified. In this 
study, we will focus on those CIRCONs that are inferred to come from AGB stars and not 
SNe. For SNe sources, it is possible to find different production zones in the models, with 
different yields, and presume that these can be sampled ad-hoc to produce the isotopic 
pattern and the carbon compounds observed. Grains that are plausibly attributed to SNe 
sources are not considered here. The problem of forming C-rich phases in standard SNe 
sources with a high O/C ratio is not really addressed. The case of AGB stars is much more 
restricted, as their envelopes are well mixed. 

As there must be a spread in the isotopic abundances in the ISM (both in time and 
space), different star formation sites will have some variability in their “bulk” isotopic 
patterns. Thus, even if we know the exact evolution of an AGB star in terms of its mass 
and metallicity, it is not possible to infer the resulting isotopic shifts. However, the isotopic 
data on CIRCONs give good insight into the nature of the nuclear effects and, often, a 
rather clear identification of the possible stellar site. We will therefore use some of the data 
on such circumstellar condensates that are associated with AGB stars in conjunction with 
processed Solar System samples to explore the implications of s-processing models of Fe 
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group elements. This requires an interplay of data from Solar System processed material 
and the data on CIRCONs. We will also use this approach to infer aspects of the isotopic 
variability in the different samples of ISM which the CIRCONs represent and to explore 
the general isotopic anomalies found in solar-system processed material. In particular, it 
will be shown that the Ti data on CIRCONs demonstrate that they show the results of 
neutron captures and that the complex patterns observed are simply a reflection of AGB 
processing on material with some small, but substantial variability in isotopic abundances 
in the original stellar sources. The questions of s-processing of Si and Ca and of variability 
within the ISM will also be discussed. 

Unlike circumstcllar condensates, the solar system was formed out of a complex 
assembly of stellar debris, from a wide variety of sources, operative at greatly different 
times, which were not completely mixed. While there is clear isotopic evidence for the 
contributions from certain precursor sources to individual macroscopic objects processed in 
the solar system, it is not possible to explain the observed isotopic variability in macroscopic 
samples by a single mechanism or a single star. What we observe are changes in the isotopic 
composition of an element relative to some terrestrial standard. This means that large 
isotopic shifts may be more easily interpretable than very small effects. The isotopic ratios 
of the “bulk” solar system is, of course, not exactly known. Positive variations in isotopic 
abundances may be easily interpreted as “additions” of an isotope to the solar reservoir; 
but the presence of negative variations clearly show that we are dealing with isotopic 
admixtures that may be different for many isotopes, as has been well recognized. 

We first present results for Cr, Fe, and Ni for AGB sources that have initial isotopic 
abundances which are solar. Then, we will explore the characteristics of s-processing 
when the source is not solar and show that there are simple rules which permit one to 
estimate relative isotopic shifts for an element independent of the mass or metallicity of 
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the particular AGB star. This leads to a general rule for stars which do not have an initial 
inventory with solar isotopic abundances. The discussion will extend to the broader issue 
of whether the observational evidence supports processing in low mass AGB stars. 


2. The Calculations 


2.1. The stellar and nucleosynthesis models 


The s-processing in AGB stars takes place on a small parcel of matter that is removed 
from the stellar envelope and is exposed to neutrons. This p rocessed materi al is then mixed 


back into the mass of the envelope present at that event flBusso et ah 


1999 : 


Kappeler et al. 


2011b As the envelope looses mass, the dilution of the processed material continuously 


changes during the thermally-pulsing (TP) phase associated with the s-process. Also the 
seed nuclei of the star are to some extent altered by the processing required to produce the 
bulk of the s-nuclei. To produce the “solar” s-inventory requires ~ 10~ 2 of all the seeds 
to be consumed. This will result in a small, but significant depletion of the seed nuclei in 
the composition of the envelope. The isotopic shifts in the envelope will be small for each 
cycle of s-processing, but they will increase in time as the envelope is greatly reduced in 
mass. While usual discussions of the s-process refer to the Fe group elements as the seeds, 
in fact there is also a major involvement of lower mass species, particularly if one considers 
elements like Ca, Ti, and Si, whose final concentration is due both to their initial inventory 
and to the modifications induced by neutron captures on lighter, abundant elements. 

The calculations presented here assumed that the initial invento ry o f all heavy elements 


was solar, with the solar isotopic ratios taken from 


Lodders et al. 


(2009), but with the 


abundances of “heavy elements” relative to H and He variable as an estim ate o f the net 


metallicity of the system (see e.g. the same procedure adopted in 


Wasserburg et al. 


2006). 














The net number of isotopes (iVj) in the envelope for all species is calculated at the end of 
each pulse-interpulse episode. This includes the s-processed material just produced and 
dredged up and the material already present in the envelope from previous episodes and not 
lost. Calculations were carried out through the full evolution of model AGB stars, including 
all stages in which any significant neutron fluence occurs. These were done up to the end of 

16 0 was followed taking 


the AGB phase. The operation of the main n eutron sourc e 13 C (« , n 
the size of the 13 C reservoir (or pocket ) from Trippella et ah (2014). A second burst of 


neutrons occurs from the 22 Ne(a,n) 25 Mg reaction in the thermal pulses; its relevance varies 
with the model, as the rate of neutron production is sensitive to temperature. We report 
on stars ranging in mass from 1.5 to 3 M & . The 3M 0 case always evolves to an envelope 
with C/O > 1, while the 1.5M 0 case becomes C-rich only at relatively low metallicity or 
after very advanced dredge-up stages. In the 3M 0 case, the neutrons from the 22 Ne source 
are important and produced at a rather high neutron density, wh ile in the lower mass 


models the 13 C source always dominates. The 13 C pocket used by 


i3r 


Trippella et al 


(2014) 


was formed adopting an exponential profile for the proton penetration down to 6 x 10 3 
M & (depth in mass). In this treatment most of the available 13 C is within about 4 x 10~ 3 


Mr..,, w hich produces the galactic enrichment of s-elcments described in 


Majorca et al. 


( 20121 ). The reactions considered in the region of Fe-peak nuclei are shown in Figured! 


where we also give approxim ate 30keV A/ 


mb, KADONIS data base 


Kappeler et al 


of interest (Takahashi & Yokoi 


axwellian-averaged values for cross sections (in 


2 0111 ) and laboratory half-lives of crucial nuclei 


1987). In the detailed calculations used, cross sections 


variations are considerable, and generally more complex than a simple 1/v scaling; also 
decay rates may change. 

The mass loss rates were taken from the new pic ture of low mass star w i nds d eriving 


from the recent works by 


Guandalini et al.l (2006) and 


Guandalini & Cristallo (12013 ). where 


the model AGB luminosities could be reconciled with observations. Cases for different 
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metallicities Z were studied; here the notations used are Z = Z 0 for a total metallicity of 


Z = 0.0154 (Lodders et ah 


20091 ): and Z = Z 0 /3. At low [Fe/H] values, all heavy elements 


are s caled downward in same way as iron. The nucleosynthesis code was that used by 

dilations presented here are somewhat different 


Trippella e t al. 


(2014); ho weve r, the ca’ 


from those reported in 


Trippella et al. 


fl20l1 as they now include explicitly 50 Cr (the full 


computer output is available from the authors on request). 


For each element, we have chosen an index nucleus (k) which is of high abundance 
and critical to any flow induced by n captures in producing isotopes of higher mass of that 
element. For chromium, iron and nickel the index isotopes were 52 Cr, 56 Fe, and 58 Ni. Note 
that some isotopes are bypassed or purely destroyed in the neutron capture chain (e.g. 54 Fe 
and 50 Cr). The production of 64 Ni depends on the complex branch at 64 Cu. The general 
reaction network is shown in Figure [Q 


2.2. Cr, Fe, and Ni in model AGB envelopes 


The enrichment factors (Ef k ) in isotopic abundances for Cr, Fe, and Ni were calculated 
for the envelopes of 1.5 and 3 M 0 model stars for Z = Z 0 and Z = Zq /3 after a different 
number of pulses; they are shown in Figure [2] and Figure [3j These calculations are for 
solar initial isotopic abundances. We define the plotted enrichment factor in the envelope 
at some particular stage (S), for the ratio of the number of nuclides i to that of the index 

nuclei k, relative to the assumed initial state (©) in the star, as: 

s 


Ef, k = 


Ni_ 

N k 


Ni_ 

N k 


o 


1 . 


( 1 ) 


We have also include d in Figures [2] an d [3] the results taken from the calculations by 


Cristallo et al. 


f201 lh. In the report by 


Cristallo et al. 


120111) the focus was on the relative 
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yield of heavy to light s-process nuclei. However, the on-line data (available from the 


http: //fruity.oa-teramo.inaf.it/) represent a 

in the region of interest ( 

Cristallo et al. 

2015) 


2015 ). They can thus be applied 


to the problem of Cr, Fe, and Ni. It is evident that both models produced marked excess 
of 54 Cr, 58 Fe, and 64 Ni. This general inference thus appears well established. There are 
clear differences found between our models and FRUITY data as can be seen by careful 
inspection of Figure [2] and [3j We note that the enrichments of 54 Cr and 58 Fe are much 
lower than that found in the FRUITY data for a 1.5M 0 source. There is larger production 
of 61 Ni in the FRUITY models in the normalized ratios shown in Figures [2] and El but the 
major excesses are always at 64 Ni. 

These issues will be noted in the text. For reasons of self consistency, the genera. 


iscussion and conclusions presented here will use the results following 


Trippella et ah 


(2014) who present higher neutron fluences and slightly higher temperatures than in the 


FRUITY models. We believe that the specific issues with these AGB models will require 
further detailed investigation. The differences found between our models and the FRUITY 
data are predominantly due to the different size of the 13 C-pocket used in the calculations 
(our pocket is about four times larger than for FRUITY), to the different solar metallicity 
assumed and to the non-linear effects induced by them. Inspection of Figures [2] and [3] 
shows in any case that the results are generally in agreement and there is no qualitative 
conflict between these independent works. In particular only the FRUITY data set and our 
results specifically treat the production of 54 Cr, 58 Fe, and 64 Ni, including the effects from 
the presence of 50 Cr and of 54 Fe, which are not relevant for computing the production of 
much heavier s-nuclei. We note that all other calculations available in the literature and 
on-line either do not have envelope abundances for Cr, Fe, and Ni, or are not adequate to 
permit inferences about yields in this atomic mass region. 
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We also notice that a full and detailed comparison with FRUITY is difficult to make, 
due to different assumptions. As an example, if we look at Ti isotopes, FRUITY outputs 
give a low production for 4 'Ti and enhancements for 49 Ti and 50 Ti at high levels (see Figure 
[4]). This is strange, considering the small pocket used in the code. H owev er, FRUITY 


adopted a total solar metallicitv of Z P: = 0.014, as in 


0.0154, adopted by us and by 


Lodders et ah 


Asplund et al. 


( 2009 ). against Z Q = 


(120091 ). This represents a 10% shift in the solar 


Z. It is possible that these small differences in the inputs cause significant variations in the 
outputs. This remains to be clarified. 


We consider that the results presented here are a reasonable effort at exploring effects 
of s-processing on the Cr, Fe, and Ni isotopes. It can be seen that the enrichment patterns 
all show clear enhancements in 54 Cr, 58 Fe, and 64 Ni and very small effects for the other Cr, 
Fe, and Ni isotopes. For later pulses the Ef 452 , Ef 8 56 , Ef 458 are quite large. It follows 
that AGB processing will produce significant shifts in the isotopic abundances of these 
three elements which resemble some of the observed patterns. The magnitude of the 
enrichment factors is dependent on the size of the 13 C-pocket, the stellar mass and the stage 
of evolution of the star. Therefore the pattern will only be approximately fixed. A more 
detailed discussion for several Fe group elements will be given later. 


2.3. General representation of isotopic results 

2.3.1. General isotope pattern from all s-process 

To gain some insight into the possible isotopic abundance patterns resulting from 
the updated neutron capture scenario outlined above, we first assume, for simplicity, that 
the composition of the envelope in a star was initially solar, and that it evolves with the 
addition of pure s-process material. The result for the envelope is simply the sum of what 
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was brought up from the neutron capture zones (n t or n k ) plus the initial inventory (TV® or 
N ’®). We then have: 


N l = N° + ni = 
n ;■ / A 7 ? 


W® 


iV® + 


T w”‘ 


Ui 


N, 


N®J fk + [n k J ^ fk) 


( 2 ) 


(3) 


where f k is the fraction of the original abundance of the isotope k remaining in the material 
processed by a pulse. The fractional shift of the isotopic ratios of the mixture relative to 
the initial state (solar, or “0”) is thus: 


ET iX = (1 - fk) 


and for two pairs of isotopes: 


T-T0, mix 
^i,k 
ttt©, mix 
h j,k 


0 ,(N?_ 

LV«J ' \N? 




/ 


A f 

N? 


- 1 



(4) 


(5) 


If Ui/rik is fixed, this implies that the ratio of the enrichment factors for two isotopes is 
constant if production ratios are constant through all stages of a stellar evolution. Note 
that fk is not the same for different elements. 

Inspection of equation (4) shows that the enrichment of a parcel of matter, once diluted 
with some original envelope material, will have an isotopic pattern that is essentially fixed 
and independent of the amount of diluting material, while the magnitude of the isotopic 
shifts will depend on fk, the fraction of the diluting material. The critical parameter 
determining the pattern can be indicated as Qf k , where: 


Q° t = 


7-1O, mix 

FI) 


riA 

n k ) 


/ 


iV® 


( 6 ) 


is the enrichment factor in the processed parcel. Insofar as Qf k is constant for all pulses, 
then this gives the isotopic abundance pattern produced by pure s-processing. The magnitude 
of the isotopic shifts in the envelope for a mixture can then be obtained as: 

ET X = QU 1 - fk) (7) 
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Table 1: Summary of the Q ?;fc values for a star of 1.5M 0 and solar metallicity, as shown by 
blue-open circles in Figure [51 Ca. Th_Cr, Fe, Ni, and Zn abundances are calculated using 


the 13 C pocket described in 


Trippella et ah 


( 2014 1 



Ca 


Ti 


Cr 


Fe 


Ni 


Zn 

i 

Qz,40 

i 

00 

i 

Qi, 52 

i 

Qi, 56 

i 

Qi, 58 

i 

Qifi 4 

40 

0.0000 

46 

3.4887 

50 

-0.565 

54 

-0.379 

58 

0.0000 

64 

0.0000 

42 

3.4598 

47 

1.3129 

52 

0.0000 

56 

0.0000 

60 

2.6721 

66 

5.9473 

43 

4.3513 

48 

0.0000 

53 

0.4613 

57 

4.9091 

61 

19.169 

67 

9.4356 

44 

1.4366 

49 

9.8874 

44 

5.4188 

58 

34.145 

62 

12.902 

68 

16.2428 

46 

8.9000 

50 

19.369 





64 

175.76 

70 

0.2361 

48 

0.2906 












To estimate Qf k we have taken, from our calculations, the outputs for the isotopic 
abundances in a parcel of matter that was removed from the envelope, subject to neutron 
exposures, following the procedure described previously, but prior to mixing back into the 
envelope. This corresponds to exploring the situation in the He-rich layers after a pulse (we 
arbitrarily chose the last pulse for a 1.5M 0 star with solar metallicity). 

The curves for Qf k for the single s-processed parcel, assuming a solar initial 


composition , are shown in Figure 0 and Tab 
a smaller 13 C pocket from the calculations of 


e [H In the figu r e we a lso show the Q values for 


Tripp ella et al. 


(1201411 (Case A) for comparison. 


The extension of this small pocket is 1 x 10 3 M 0 and it contains 3.1 x 10 6 M 0 of 13 C. 


We note that the Q values for 54 Cr, 58 Fe, and 64 Ni are large as compared to other 
isotopes of the same elements. They are therefore not sensitive to the initial condition. The 
Q values are directly dependent on the size of the 13 C pocket, so that this parameter must 
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be taken in full consideration. If we restrict our discussion to a fixed pocket size, then some 
general rules apply. General patterns will be approximately the same for different pocket 
sizes, but will be shifted upwards or downwards as can be seen in Figure 0 (lower curves). 


We now compare the results for Q t , k to the pattern obtained in the full calculations of 


E l k in the envelope, for our standard 13 C pocket as given i 


m 


Trippella et al. 


(120141 1. These 


enrichment factors are shown in Figure [6] and may be compared with the approximation 
represented by Qi )k shown in Figure [5j All the results are normalized to the isotope with 
the largest shift for each element at different stages. This exhibits a comparison of the 
isotopic abundance pattern for each element as a result of the full calculation. The results 
for other Fe group elements are also shown in Figure [6] It is clear that all the s-processing, 
almost independent of stage or metallicity, gives essentially the same isotopic pattern, but 
with different magnitudes of the enrichment for a fixed 13 C pocket. Further, this pattern is 
almost indistinguishable from the quantitative pattern approximately obtained by using the 
Qf k reference values. We can now explore the dependence of on the role of the neutron 
fluence. 


We used a reference state for Q for the last pulse of a 1.5M 0 , Z = Z Q model AGB. 
First, note that a key term in the denominator in Q i)k is n k , the number of index nuclei 
returned to the envelope by the processed parcel. For Fe, the index isotope is also effectively 
the remaining 56 Fe inventory in the processed parcel. As 56 Fe is not produced from lower 
mass nuclei, it is only destroyed. Higher neutron fluences cause a decrease in 56 Fe in the 
returning parcel. Hence, an increase in the fluence would mean a decrease in Qj 56 . For Ni, 
the index isotope is 58 Ni: it is mainly destroyed, but also produced by inflow from neutron 
captures on Fe (which produce more Ni). Thus, there is a much enhanced increase in 64 Ni 
as compared with 58 Fe. Similar arguments follow for Ti and_Cr. In Figure [5] we compare 


our results for a large 13 C pocket (like that used by 


Trippella et al. 


2014) with those from a 
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smal 


2006 


er 13 C reservoir of 10 3 M 0 (as discussed e.g. in lGallino et al.lll998l: IWasserburg et al. 


The effects of s-processing on the parcel considered clearly define the isotopic patterns 
for all models explored. It appears that the above approximation is robust and reasonably 
accurate. It is importance to recognize that details for each pulse are not required to 
describe the resulting isotopic enrichment pattern. 

It follows that, for a given 13 C pocket, the s-process produces relative isotopic shifts of 
an element which are, to a good approximation, independent of the metallicity, model or 
stage of the particular AGB star, provided the star had initially solar isotopic ratios for the 
elements of interest. 


2.3.2. The Relative Magnitudes of the Enrichment Factors E^j 

While the isotopic abundance pattern in the envelope is almost an invariant for all s 
processing models, the magnitude of the effects for each element is very sensitive to the 
number of pulses, the efficiency of TDU,the abundance of 14 N in the processed parcel, 
metallicity and the stellar mass. It is also dependent on the size of the 13 C-pocket used. 

In Table 2 we show the relative magnitudes of E i)3 at the last pulse for Ti, Cr, Fe, and Ni 
in models of different mass and metallicity, but with Solar initial isotopic ratios. Only the 
isotopes with large shifts in patterns are considered here ( 50 Ti, 58 Fe, 64 Ni). Data are taken 
from the inserts in Figure [6l 

It can be seen that the relative effects for these isotopes with major shifts in patterns 
relative to 50 Ti are very different. 54 Cr is enhanced much less than 50 Ti in all cases. 
However, 58 Fe is enhanced as much as or more than 50 Ti. For 64 Ni, the enhancement is 
much greater than for 50 Ti. It appears that this is a general rule. The relative effects to be 
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observed depend on the scaling between different elements. In consideration of variability 
in the initial isotopic composition, all samples from AGB sources with few pulses will not 
show significant or recognizable effects. Stars formed at early epochs from low metallicity 
clouds will show the largest effects. 


2.3.3. Measured Isotopic Shifts and the Initial Composition of a Stellar Source 


It is widely recognized that the initial isotopic compositions for stellar sources which 
formed from the ISM are variable. To discuss the isotopic effects resulting from AGB 
processing of initially non-solar isotopic patterns, we must explore the consequences of 
such variability. The numerical calculations typically assume solar isotopic ratios. In 
consideration of variability in the initial isotopic composition, all samples from AGB sources 
with few pulses will not show signicant or recognizable effects. Only those samples of stars 
with advanced pulses will show clear signs of s-processing. Low metallicity from earlier 
events will show the largest shifts. 


In case of a non-solar initial composition, we define the enrichment factor in the 
envelope at some particular stage (A), for the ratio of the number of nuclides i to that of 
the index nuclei k, relative to an assumed initial ( non-solar ) state “O” in the star, as: 


E? k = 


N k 


N k 


o 


- 1 


( 8 ) 


Again, when Ef k < 0, this corresponds to a depletion of the isotope i relative to k. 
If the initial ratios in the star are (Ni/N k )° , then we define 




(9) 


As we assume that Ui/n k in the returning processed parcel is the same for all s processing 
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in an AGB star, the enrichment rule stays the same. Hence, equations from (2) to (5) can 
be generalized to an arbitrary initial state (0): 

E°„ = (1 - ft) [(n,/nO / (N°/N°) - l] = (1 - ft) Q° k (10) 

Further, as rii/n k is taken as fixed, then: 

Q?,k (1 + a i,k) = Qf )k — ®i,k (11) 

When measurements are made of stellar debris, the standard reference state for isotopic 
ratios is solar. Thus all measurements so represented are () 0 , independent 
of the actual initial value in the source. The effective calculated enrichment factors are thus: 


E% lc = (1 + Ci,t) E°r x + a 


i,k 


( 12 ) 


It follows that: 


( TT'Calc 
\ L i,k 


“«) / ( E f,t c ~ “**) = {Q?,t ~ “a-) / (Qf.t - a i.t) 


(13) 


This is the explicit, general relationship between the calculated enrichments from laboratory 
measurements to any possible AGB source, assuming the approximations to be adequate. 

It relates the measured isotopic ratios to a standard reference s process source Qf k , as given 
above and to the shifts in the initial composition of the source. Note that there are as many 
values of as there are isotopic pairs, so the system is not determinate. However, if one 
of the ck’s can be estimated, then all of the others may be computed from the calculated 
delta values for a sample and some standard Qf k , such as given in Table [0 


Comparison of stellar models with results for a sample (A") measured in the laboratory 
requires some special care. The range of effects is enormous and the measured laboratory 
results are presented in different forms. They are given in units of 5 (parts per mil) or e 
(parts per 10 4 ) as defined below: 


= 10 a ( /; 


J i,k 


j-Sample / t-q 


l i,k 


AS 


— i n-1 


4 


(14) 
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so that: 


E O.Cic = w - H X k = 10 -4 £ X 


(15) 


Equation (13) can be then rewritten, in terms of “delta values”, as: 


(lO- 3 Ct iC - “<•*) / (lO- 3 -*" 0 - “m) = (Q?,t - «.,*) / (Q% - a**) 


(16) 


Here is the ratio of the ion beams of isotopes i to k measured in a sample, corrected 
by a discrimination factor (1+Zb;^). For ion-probe data, the discrimination factor is 
typically determined from the measured ion beam ratio lf k in a terrestrial sample. It is 
usual, in many cases, that instrumental discrimination factors are determined by one pair 
of isotopes (w, v ) relative to the nominal Solar standard (usually a terrestrial material) 
and appropriately applied to other isotope pairs assuming some functional dependence on 
mass. It is important to note that the experimental errors for the delta representation 
used for Sputtered Ion Mass Spectrometer (SIMS) or nano-scale SIMS (nano-SIMS) are 
typically given as one sigma. In contrast, the smaller effects measured by Thermal Ion 
Mass Spectrometry (TIMS) or Ion Coupled Plasma Mass Spectrometry (ICPMS) are two 
sigma errors. These effects are much smaller, but are more precise. They often represent 
a normalization to a terrestrial ratio for two isotopes to correct for in-run fractionation. 
When isotopic ratios are shifted, the normalization can lead to spurious effects. 


The problem to be addressed is to relate a set of observations for many elements and 
isotopes in an individual CIRCON to the composition of a single stellar envelope at some 


Z inner 


2014). That is, one seeks 


stage “S” that would produce such a grain (see review by 

G) S ‘ 

to find a stellar model where at some stage Ef k is related to the observed value in the 
grain. 


In this case of laboratory measurements of an exotic material, we take Eff alc from 
equation (15). For the case of processed, but not completely homogenized samples of solar 
system material, one tries to relate the effective average of some sources (Ef k ) to the 
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observed isotopic composition by some undetermined scaling factor and seek to infer which 
sources are plausible. Large effects for different elements may be more definitive on a single 
object. In all cases, the variation in the initial state of stars relative to the Solar reference 
values must be considered. 

We will first, for simplicity, calculate the nature and magnitude of effects for Cr, Fe, 
and Ni and then discuss the other Fe group elements. We will initially assume in our 
numerical calculations: /N° = Nf /N® (which is not in general true). 

We will show the high degree of enrichment for 54 Cr, 58 Fe, and 64 Ni. For Ti, a deformed 
J-shaped trend is exhibited, with a major effect in 50 Ti. 

As the numerical calculations used solar values, we must consider the extent to which 
variations in the original stellar source would affect the results. From Figure 0 we note 
that the Q lyJ values for 54 Cr, 58 Fe, and 64 Ni are large, but are small for all the other Cr, 
Fe, Ni isotopes. As a result, changes in due to changes in the initial composition (say, 
a^k ~ 2) will not greatly alter the general isotopic pattern, with the dominant effects being 
at 54 Cr, 58 Fe, and 64 Ni. Hence, the resulting patterns for Cr, Fe, and Ni are robust and will 
be roughly fixed for the AGB models considered here, independent of small shifts in the 
stellar initial abundances. 


3. Full isotopic results for AGB stars: Cr, Fe, and Ni products 

In this section we present the detailed results for the isotopic composition of the 
envelope, represented by Ef k , over all pertinent stages of a model star evolution, from the 
first formation of a 13 C pocket until the last s-process event, prior to the transition of the 
star into a planetary nebula (see Figure [6]), which leaves a white-dwarf remnant. These 
are the results of the full stellar model calculation and not the partial ones, referring to 
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the He shell, given above in Figure [5j It is readily seen that the patterns for all cases are 
very close to that obtained for a single s-processed parcel (see Figure El), but that there 
are small to significant shifts depending on the details of the particular stellar mass and 
dredge-up episode. We will discuss the effects of different initial states later in section 5. 
The degree of enhancement increases with the number of pulses due the repeated neutron 
processing and also to the decrease in the remaining mass of the envelope. In Figure El 
we plot the data normalized so that E % ^ is set to one for a relevant normalization isotope 
(e.g. the isotope 64 Ni for nickel), using for reference an isotope that is much enhanced. 
The figure inset gives the scaling factors for each case. Note that there are a wide range of 
scaling factors, but that in general, the forms of the curves for all cases are similar. This 
is a reflection of the nature of s-processing in the Fe-peak, as shown in section 2.3.1, and 
which generates isotopic shifts for Cr, Fe, and Ni with clear major excesses in 54 Cr, 58 Fe, 
and 64 Ni relative to the index isotopes and only small effects for the other isotopes. The 
circumstcllar condensates with Cr, Fe, and Ni formed around AGB stars will have a wide 
range of enrichments (as revealed by the scaling factors in Figure [6]) depending on the stage 
of mass loss and dust formation and on the stellar mass and metallicity. 

Note in particular that an enrichment of 54 Cr/ 52 Cr is found that can run up to 
-^ 54,52 ~ 1 for the bulk envelope at late stages, for the low metallicity case. We notice 
that, if some material is ejected without full dilution in the envelope, at any stage, then 
the effects can be very much larger, with the isotopic pattern being fixed by the purely 
s-processed parcels discussed in the previous section. 


4. Meteoritic data and comparison with model predictions 

It is important to note that, while excesses in an isotopic ratio relative to the solar ratio 
may be explained by addition of that isotope to a Solar mix, negative values require that 
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Table 2: Ratios of E itk for 54 Cr, 58 Fe, 64 Ni, 42 Ca, and 29 Si with respect to 50 Ti. Calculations 
refer to the last thermal-pulse of 1.5 and 3M 0 stars assuming two different metallicities (solar 
and 1/3Z 0 ). 


Ratio Element Ratio 

1.5 M e 

Z Q Z & /3 

3 M q 

Zq Zq/3 

-^(54,52) /E(50, 48) 

Cr/Ti 

0.35 

0.16 

0.49 

0.27 

^(58,56)/A'(50,48) 

Fe/Ti 

1.65 

1.00 

3.57 

1.90 

-^(64,58) / -^(50,48) 

Ni/Ti 

7.49 

2.16 

5.33 

3.08 

-^(42,40) /^(50,48) 

Ca/Ti 

0.22 

0.14 

0.22 

0.18 

^(29,28)/^(50,48) 

Si/Ti 

0.14 

0.05 

0.09 

0.03 


something is being added having an isotopic abundance pattern for several isotopes which 
is quite distinct from the Solar pattern. This is a well known fact that is often ignored. 
For all ion-probe data on CIRCONs where the isotopic shifts are large (many hundreds 
per mil), we believe that all the data processing that has been applied is to correct for 
instrumental discrimination effects, by using the isotopic ratios of an isotope i relative to 
some index isotope k through the determination of that ratio in a terrestrial standard for all 
isotope pairs. This correction is then applied to the measurements of all data. In contrast, 
the data on CAIs, which are far more precise (e units), are most often normalized to a pair 
of isotopes to correct for instrumental fractionation. This procedure effectively requires one 
pair of isotopes to have the terrestrial ratio and thus affects the pattern in all the other 
ratios. This difference in the treatment of the experimental data must be kept in mind. 



- 22 - 

4.1. Cr 


The discovery of i sotopic anomalies in Cr by 


Birck fe Allegro (! 198 1) and 


Birck fe bugmairl (1988) established the presence of small 54 Cr/ 52 Cr ex cesses (of the 
order of ~ 6 x 10 -4 , or 6 e units) in CAIs from the Allende meteorite. Pap anastass iou 


(119861 ) then measured Cr in the two “fractionation and unknown nuclear anomalies” (the 


so-called FUN) inclusions, (EK1-4-1 and C-l) that had previously been discovered to have 
very large isotopic anomalies in many elements. All these mass spectrometric data are 
normalized to 50 Cr/ 52 Cr to match a terrestrial standard. Considering the magnitude of 
isotopic shifts in CAIs, he found a very large excess in 54 Cr (48 e units) and a large excess 
in 53 Cr (16 e units) in EK1-4-1. However, for C-l the effect is a large depletion in 54 Cr (—24 
e units) and data on non-FUN samples displayed subdued effects that were intermediate in 
behaviour (non-FUN samples of CAIs are also shown in Figure 1 of that paper). 


The issue is then to understand whether any of these patterns can be plausibly 
explained by an AGB model. In general, most CAIs exhibit excesses in 54 Cr. Analyses 


of bulk meteorite samples (CV, C and CM types ' 


level. However, an investigation by 


Rotaru et al. 


exhibit no effects at about 1 e unit 


19921 ) showed that these meteorites, if 


subjected to acid leaching, yielded widely variable Cr isotopic compositions in any given 
leach sequence on a single meteorite. They found, for the Cl chondrite Orgueil, a range 

eaehes. In most cases, the results 


of values from ~ —9 to ~ +100 for e m in different 

showed strong excesses of 54 Cr. Pod os ek et al.l (1997) sought to find the carrier grains of 


anomalous Cr by a stepwise dissolution process and found 54 Cr excesses of up to ~ 210 e 
units, but no identification of carrier grains. It is clear that there are materials of highly 
divergent nuclear sources in any bulk sample^ Several groups have carried out searches for 


the carrier grains of the anomalous Cr. 


Dauphas et al. 


(120101 ) presented an extensive study 


searching disaggregated and treated material for carrier grains from 1,400 to 30 nm in size. 
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There were clear excess in ° 4 Cr, but not of very large amounts (up to 17 e units). Only one 


exceptional grain with E 5 ^ 52 = 2.9 ± 0.33 was found, which was attr ibutec 


source on account of its very large excess. In a study by 


Qin et al.l (120111) on insoluble 


to a supernova 


residues with Cr-rich grains (~ 200 nm), they found three grains that are enriched up to 
1.5 times the solar value, and concluded that this must be the result of supernova debris 
injected into the protosolar cloud. These grains are plausibly circumstellar condensates; 
however, there are no measurements made on identified circumstellar grains that are clearly 
related to AGB sources; the measurements themselves are very difficult to perform and 
have large errors. The Cr data obviously show clear effects, with frequent 54 Cr excesses and 
some 54 Cr deficits; but full measurements for all Cr isotopes are limited and no data are 
available on well-defined circumstellar condensates that are clearly related to AGB stars. 


4.2. Fe 


Extensive data on a variety of bulk meteorites and some Fe-Mg chondrules from 
Allende show n o shifts in the Fe isotopes, inc l uding 58 Fe. The only available data on Fe 


anomalies is by 


solutions of EK1-141 used for Cr measurements (Pap anastass i ou 


Voelken ing & Pa panastassiou (1990). The results^ obtained on the same 


1986 ) are of high precision 


and normalized to the solar value of 54 Fe/ 56 Fe. They show clear, large excesses of 58 Fe of 
292 e units and no shifts for 57 Fe. The C-l sample shows instead only marginal enrichments 
in 58 Fe; other CAIs exhibit small, but cle ar excesse s in 58 Fe. The data on Fe in CIRCONs 


are limited. The first reports were by 


Hoppe et ah 


was on AT-grains that are associated with SNe. 


( 2000 ). 


Hoppe et aL- 


Using an IMS-3F, their focus 


2000) only measured 54 Fe/ 56 Fe 


and did not detect any effects. These are difficult measurements and_they noted 


to several isobaric interferences including 54 Cr. A recent report by 


Marhas et al 


imits due 


( 120081 ) 


presented extensive data on Fe and Ni isotopes in SiC grains using a Carneca NanoSims. 
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Again, there are limits on which isotopes conlcl be measured dne to isobaric interferences, 
in particular for 54 Cr, 58 Fe, and 64 Ni. As a result only 54 Fe/ 56 Fe, 57 Fe/ 56 Fe, 60 Ni/ 58 Ni, 
61 Ni/ 58 Ni, and 62 Ni/ 58 Ni were measured with corrections for 54 Cr to 54 Fe and for 58 Fe to 
58 Ni, assuming terrestrial values for the corrections. The focus of this work was again on 
X-grains associated with SNe. Large excesses were found for 57 Fe and 61 Ni, 62 Ni in some 
X grains. Extensive measurements were also made on Mainstream grains associated with 
AGB stars. These showed no large shifts (see their Figures 4 and 5). The major emphasis 


m 


Marhas et ah 


(120081 ) was to identify the effects they measured with selected zones in a 


reference SNell model. This approach assumes that different zones are sampled to produce 
the isotopic pattern and impart it in to a SiC grain. These workers also presented the 
effects of s-processing Fe and Ni in the envelopes of AGB stars. However, they present no 
calculations for 54 Cr, 58 Fe or 64 Ni. We presume that these were ignored because they could 
not be measured. As we showed in section 2, these are the isotopes that will exhibit, by 
far, the largest effects from s-pr ocess ing in AGB stars; as a consequence, the results on 


main stream grains in the 


Marhas et al. 


(2008) report do not either support or negate the 


findings presented in the present work. 


4.3. Ni 


The discovery of small (a few e units) anomalies in Ni was found bvlBi rck fe Lugma ir 


( 19881) . The l argest e ffect was on 64 Ni (about 3 e units). Measurements of bulk meteorites 


(IRegelous et ah 


20081 ) showed no isotopic shifts at a high level of precision (< 0.1 e units) 


for many meteorites, but no data on 64 Ni were present. 


An important, critical measurement on Ni was made by 


Chen & Pa panastassiou i 


(2014). These workers showed endemic 64 Ni effects in CAIs which ranged from +55 e units 
(Egg-3) to very low levels (~ 1 e units) in some other samples. The data were normalized 
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to the 58 Ni/ 61 Ni value for terrestrial Ni. The effects on all other isotopes were small. One 
sample showed a distinctly negative value (~ 1 e) at 60 Ni (with their normalization). The 
observation on Egg-3 was the focus of our interest as it showed unambiguous nuclear effects 


with a clear pattern for this key “Fe group” elemenh Again,_t 


on CIRCONs, other than the results reported by 


le re ar e no Ni isotopic data 


Marhas et al 


( 2008! ). which, as indicated 


above, are not adequate for the problem under discussion. 


4.4. Ti 


There are abundant isotopic data on Ti, for solar-system processed samples (CAIs and 
many bulk meteorites). The patterns are complex and varied. Many publications on Ti 
isotopes in processed solar system material are available, but the results we mention below 
are typical of the whole range observed. It is often found that there are well defined excesses 
at 49 Ti, 4 'Ti, and 50 Ti, with the 50 Ti anomalies being larger than for 49 Ti (as found in 
EK-1-4-1). There are also samples with distinct, large deficits in 50 Ti and smaller, but clear, 


deficits in 49 Ti (C-l): (see 


Niederer et ah 


198(1 with a normalization to solar 48 Tl/ 46 Ti). 


Absolute isotopic ratios were obtained using a, modified tec’ 


normalization for instrumental fractionation (Ni edere r et al. 


mique that eliminates the 


19851 ). These workers showed 


that at least three, and possibly four stellar sources must have contributed to the five Ti 
isotopes to make up the samples of processed materials. An extensive discussion of possible 
stellar sites and sources was also presented by them in trying to identify the compositions 
of the exotic components. It is evident that the processed solar materials do not yield a 
direct tie in with individual possible stellar sources. 


Because of the problem of complex mixed sources, we turn to the Ti data in CIRCONs, 
with the hope of finding some regularity using known nucleosynthetic signatures of AGB 
stars. The AGB processing for a single grain then depends on the stellar evolution and 
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the initial composition of the parent star. In partic ular, the variations observed in Ti 


isotopes found in oxygen-rich CIRCONs (Choi et al 


1998) suggested that the Ti isotopic 


compositions follows the pa' 
by Gallino and co-workers. 


tern gen erate d by AGB processing, following on the calculations 


Nittlcr et al. 


(120081) presented an extensive study of Al, Ca, and 


Ti in CIRCONs. They showed very strong evidence that most of these grains were formed 
in red giant and AGB parent stars. They reported Ti data on three oxide grains that 
showed small isotopic effects, but with divergent patterns. These grains must have been 
produced in a regime where C/0<1, implying a low mass star prior to C enhancement. 


Extensiv e studie s of the isotopic composition of many elements in SiC circum s 
^OlJ ) were presente d by several workers ; see e.g. 


cellar 


condensates (Zinner 


in particular, see their Figure 5. See also 


Amari et al. 


Amari et ah 


(12001D; 


(120041 ). in particular see their Figure 


4. The data exhibit a wide and diverse array of very different patterns; these last have 
often been described as being “V”-shaped, but with large variations, sometimes negative. 
We will rather use a “J” shape as a descriptor, because that is the pattern which would 
be produced in AGBs from an initially Solar mix (see Figure [7j, in seeking some hint as 


to which examples from the zoo of Ti isotopic patterns migh 


that several samples show “J” shape patterns in 


Amari et al. 


serve as a guide. We note 


(120011) ; the same is true for 


Huss fe Sm ith (120071 ). see in particular their Figures 2, 3, and 11. The patterns of many SiC 


grains are quantitatively indistinguishable from those of AGB models for an initial solar 


Amari et al 


020011 ) for KJGM1-158-5, 


isotopic abundance. We show the measurements by 
which has a very large enrichment in 50 Ti, in Figure 0 For comparison, we include the 
results for Q it48 solar, normalized to the measured A 50 ,48 value. It can be seen that within 
the errors, the results are indistinguishable, if one considers the errors in both 50 Ti/ 48 Ti 
and 49 Ti/ 48 Ti. 


We now consider a wider number of samples, some of which have only more modest 
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effects in £ 50 ,48 (down to 200 per mil). 


101, 58b, and 110 from 


Using data from 


Arnari et al. 


(2001 


2004) and grains 


Huss & Smith (2007), we show in Figure El the results, normalized 


to the measured £ 50 , 48 , including also error bars. Again a set of normalized reference curves 
using the Qi^s values are shown. Some samples, such as KJGM4C-314-3, appear to be 
widely deviant with a 4!, Ti shift above that of 50 Ti. However, if one considers the one-sigma 
error bars of 20% for £ 49,46 and of 24% for £ 50 , 46 , there is no clear discrepancy. From 
these observations we conclude that a large fraction of these grains show clear evidence 
of s-processing. However, there are many cases that also show clear disagreement with 
significant s-processing. These differences, shown in Figure El would correspond to shifts in 
initial isotopic abundances of 100 per mil for some isotopes which are not enriched in the 
s-process calculations presented here. 


We are then left with the issue of what the variations are due to. Many of the 
samples showing clear deviations from AGB models still exhibit aspects of the “J” or 
“check” pattern. Note that for 4(l Ti (which is 8% of the total Ti), a 5% increase in the 
46 Ti/ 48 Ti ratio would give (relative to Solar) a 5% shift in 46 Ti/ 48 Ti and change the 
“J” pattern produced by s-processing into a “V” pattern. There is a very strong indication 
that many of the variations observed represent s-processing, but from somewhat different 
initial isotopic compositions. When the s-process effects are large (due to grain formation 
in very advanced evolutionary stages), the Ti isotopic pattern closely approaches that of 
pure s-processing. When instead grains form at earlier dredge-up episodes, then the initial 
composition dominates. 


4.5. Si 


We now consider Si isotopic shifts in CIRCONs in the light of the above discussion 
of Ti which shows clear evidence of s-processing, but with some intrinsic variability in 
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the initial values due to granularity in the source ISM. The data on SiC grains (main 
stream and A, B grains) associated with AGB stars exhibi t a c l ear t rend in the isotopic 


Zinner 


2014). The fundamental 


ratios of Si, as found by other workers (see Figure 3 in 
iss ues o n the isotopic composition of Si in the ISM were fully recognized and discussed 


by 


Timmes & Clayton (1996). These workers discussed the diversity of SNe sources and 


the possibility that these effects could be used to trace earlier galactic evolution. The 
observed trend expressed in terms of 5f 0 28 versus (*> 29,28 ^ ias a sl°P e close to (or larger than) 
unity, with some points displaced in the direction of a line of much smaller slope (~ 0.3). 
The slope of ~ 0.3 is calculated for a model with 1/v dependence for neutron capture 
cross sections. This inconsistency of predictions with the main trend of high slope was 
recognized by many authors. Some significant efforts have been made to explain these 
results. We note that the cross sections do not follow the 1/v rule and are subject to large 
uncertainties, especially for 30 Si (Franz Kappeler & Iris Dillman, private communication). 
Hence the slope predicted for AGB stars may be more complex than previously assumed. It 
is dependent on the effectiveness of 13 C burning and on the extension of the 13 C reservoir. 
For small pockets, the slope is small (< 0.5) and determined mainly by neutrons from the 
22 Ne(a,n) 25 Mg reaction. For large extensions of the pocket the slope increases substantially, 
up to about 0.8 in our models. In no case can large enrichment factors for 29 Si and 30 Si be 
produced in AGB envelopes, as the 28 Si neutron capture cross section is small (this is also 
the case for Ca isotopes). While more advanced AGB computations may change the results, 
we consider here that the SiC measurements are plausibly explained by other important 
factors. In particular, the possible intrinsic local variability in Si in the local ISM from 
which the grains formed must be considered. It is of importance to note that the relative 
enrichment of Si and Ca effects from AGB processing as compared with Ti, here is rather 
low (see Table [2]). In contrast, 58 Fe and 64 Ni should exhibit larger effects than 50 Ti. For 54 Cr 
there should be smaller effects than for 50 Ti. Thus, in consideration of the small variability 
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(10% to 20%) in the ISM isotopic composition estimated from the Ti data, we expect that 
the coupling of Si data with Ti due to AGB processing will be very weak. The dominant 
effect for Si is due to intrinsic variability in the contributions from distinctive SNe sources. 
Coupling the results for both Si and Ti caused many workers to invoke a model of Galactic 


Chemical Evolution to interpret the observations. The original proposal 
isotopic data in CIRCONs by Galactic Chemical Evolution was made by 


;o explain the O 


Nittler fc Coswik 


(1 99 7) and discussed further by Nittler (2005) in a stochastic model, considering mixing 
from different sources. 

Efforts to relate the Ti and Si effects to Galactic Chemical Evolution resulted in some 
difficulties. As pointed out by 


Amari et al. 


(I200lh and by [Nittl er (20051), the correlation 


lines for Si and Ti isotopes for the Galactic Evolution trend do not pass through the solar 


value. Further, Nittler (l2005h pointed out that the Si trend could be explained if the role 
of SNela ejecta in p roducing 28 Si was la rger than assumed in the previous years. In an 


important paper by 


Lugaro et al. 


(119991 ). after the work of 


Timmes & Clayton (1996), the 


following points were made: i) that there was local heterogeneity in the ISM at the time 
when the “main stream” grains were formed; and ii) that there are shifts to be expected in 
the Si isotopic abundances in the ISM due to intrinsic variability in the ratios produced by 
core collapse SNe and by SNela. These workers still argued that there are effects of Galactic 
Chemical Evolution, but considered the possibility that there is intrinsic heterogeneity in 
the ISM at essentially all times. 


The wide possible range in 29 Si/ 28 Si and 30 Si/ 28 Si presen t ed by 


was based on the detailed SNe yielc 


subsequent report by 


Rauscher et al. 


s by 


Woosl ev & Wea ver (119951 1 and is seen in a 


Lugaro et al. 


(1999) 


(12002 1. We note that recently published yields of the 
Si isotopes for SNela clearly indic ate that, indeed, 29 Si and 30 Si ar e grossly under produced 


as compared to " 8 Si in the ejecta (Bravo & Martinez-Pinedo 


2012 ). SNela are major source 











































30 


of Fe. The above authors showed that 29 Si/ 28 Si ~ 1.96xl0~ 3 and 30 Si/ 28 Si ~ 3.6xl0 -3 . 
These are about 1 / 10th the solar values. The slope of a line connecting Solar and these 
SNela values is 0.93 in the 5f 028 versus $^ 2 s- 

About 15—45% of the 28 Si in the bulk Solar inventory is from SNela sources 


(JLugaro et al. 


19991) and all of the remainder is from SNell. The various SNell sou rces 


(depending on mass) show a wide spread in Si isotopic abundances flWooslev fe Weaver 


1995 


Rauscher et al 


20021) . Calculating a grand canonical average using an initial mass 


function that gives the large scale Solar Si isotopic ratios provides a plausible base line as 
a guide. However, assuming Solar abundances as a universal rule is not plausible. It is 
evident that the bulk Solar value is the result of a blend of di 
29 Si and 30 Si coming from very massive stars (> 30M 0 , see 


iverse sources wit 

l much of the 

Chieffi & Limongi 

2013 

). If we 


consider inhomogeneities in the ISM resulting from rather small variations in the mixing 
ratios, then a 1-10% change in the contribution of SNela sources to a local star formation 


region would result in a shift of 100 per mil for Si isotopes. This is esse ntially the 
that is observed for the mainstream and the A and B grains (see 


is essentially 

the 

Amari et al. 

2001 


rend 

and 


references therein). However, one must further note that the Si data show that the Y grains 
and the Z grains are displaced to the right of the main array (which has a slope of ~ 1). It 
is possible that this displacement is due to AGB processing of materials that were on the 
slope ~ 1 line. In general, the X grains (which have very negative deltas) are recognized as 
having rather direct connections with SNe condensates. 


The aforementioned study by 


Nittlerl (2005) of a stochastic sampling of SNe ejecta 


from different sources found that this approach does not reproduce the strong correlation 
between 29 Si/ 28 Si and 46 Ti/ 48 Ti ratios observed in CIRCON grains. However, this approach 
attributes both Si and Ti isotopic effects to SNe. If the dominant effects for Ti in 
mainstream CIRCONs are due to AGB re-processing and the Si effects are from SNe, then 
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this correlation is not applicable. The O isotopic patterns are in many cases not well 
understood and may indeed require SNe contributions. More work is needed on this subject 
with full attention given to intrinsic heterogeneities in the ISM. 

We infer that the Si isotopic data of CIRCONs found in meteorites are plausibly 
explained by modest variability in the isotopic composition of Si in star forming regions. 


4.6. Ca 


Measurements of all the Ca isotopes are quite difficult (particularly for the rare isotope 
46 Ca, which is only about 3 x 10~ 2 of the total Ca). There are also isobaric interferences 
from 46 Ti and 48 Ti for the corresponding mass numbers. Non-linear isotopic anomalies of 


Ca (n ot attr ibutab 


by 


Lee et ah 


e to mass dependent isotopic fractionation) in CAIs were first discovered 


( 1978 ). Normalizing their data to solar 40 Ca/ 44 Ca, they found that 48 Ca was 


enriched in EK141 by +13.6 per mil. For sa mple C-l they found a defici t of ~ —3.0 per mil 


in 48 Ca. No effects were detectable in 46 Ca. 


Ni edere r fe Papanastassio ul (1198411 determined 


absolute isotopic ratios for CAIs, eliminating the normalization procedure to correct for 
instrumental fractionation effects. They confirmed the effects reported for EK141 and Cl; 
they also fo und a CAI with very large excesses in 46 Ca (~ 7 per mil) and with no effects 


m 


42 ’ 43 ’ 48 Ca. 


Papanastassiou fe Brighaml (19891), in a study of Mg, Ti, Cr, and Ca in CAIs, 
found clear, very large deficits in 48 Ca (~ 31 per mil) and no detectable effects in 46 Ca. 
These samples also showed large correlated deficits in 19 Ti, 50 Ti, and 53 Cr (see their Figure 
1). These reports clearly demonstrate significant isotopic effects in Ca in processed solar 
materials and represent the l argest eff ects so far found in them. A recent work of very high 


precision was carried out by 


Huan g et ah 


(2012) on CAIs and found that there were no 


detectable shifts in 43 Ca/ 44 Ca or 40 Ca/ 44 Ca (below ~ 1 e units), but several samples showed 
clear excesses of 48 Ca (~ several e units). In their study, these authors used a double 
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spike technique so that some of the data are not normalized to correct for instrumental 
fractionation. The absolute 44 Ca/ 40 Ca and 44 Ca/ 42 Ca data are very precise (±0.06 per 
mil). However, the 48 Ca/ 44 Ca ratios were calculated assuming the terrestrial value. This 
means that even the small errors in 44 Ca/ 40 Ca and 44 Ca/ 42 Ca would propagate into shifts 
of over 1.2 e units in the 48 Ca/ 44 Ca ratio, which is not much less than the effects reported 
considering the listed errors. No measurements were made of 46 Ca for the reasons given 
above. They also found clear evi dence of mass-dependent f ractionation in both C AIs and 


bulk meteorite samples (see also 


Huang & Jacobsen 


2012 ). 


Movnier et ah 


(.2010) leached 


samples of a C-l chondrite searching for carrier grains of 54 Cr and also measured Ca. They 
found widespread small (~ 1 e units) variations in 40 Ca and examples of 46 Ca deficits in 
one case and a couple of cases where there was clear small 48 Ca excess (~ 3 e units). It is 
important to note that large errors are often found for 46 Ca measurements^ They further 


H-W. Chen et ah 


noted that 54 Cr effects were not correlated with those on 48 Ca. 
carried out an extensive study of 48 Ca anomalies in a wide variety of bulk c 


( 2011 ) 


meteorites, again using a normalization to 42 Ca/ 44 Ca following 


Russell et ah 


ifferentiated 


119781). They 


confirmed the ubiquitous, but small, isotopic heterogeneity in planetary bodies of the Solar 
system for 48 Ca and also reported one sample with 46 Ca excess far outside their typical 
experimental error. In their discussion, they focused on SNe sources for the 48 Ca effects, 
since 48 Ca is produced in SNe. 


The isotopic data on Ca for CIRCONs are quite limited. There are no measurements of 
46 Ca and very few for 48 Ca. Note that sputtering ion mass spectrometry, which does not use 
preparatory chemistry to remove interfering elements, requires that interferences must be 
resolvable and identifiable. The interferences from 46 Ti and 48 Ti make such measurements 
very difficult for a rare isotope like 46 Ca. 

measured 42 Ca/ 40 Ca and 44 Ca/ 40 Ca in two CIRCON oxide grains 


Choi et ah 


119981 
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and found them to be the same as the solar value, but with large uncertainties. 


Choi et al 


(1999) subsequently measured two more CIRCON oxide grains and found 42 Ca/ 40 Ca, 


43 r 


Ca/ 40 Ca, and 44 Ca/ 40 Ca to be normal within errors. 


Nittler et al. 


(2008) presented a 


very extensive and thoughtful study of Al, Ca, and Ti in a large number of CIRCON oxide 
grains associated with AGB sources. They presented data on 42 Ca/ 40 Ca, 43 Ca/ 40 Ca, and 
44 Ca/ 40 Ca. In five cases they found shifts in 42 Ca/ 40 Ca from +249 to —38 per mil relative 
to solar. For 43 Ca/ 40 Ca they found shifts of +118 per mil in one sample and for 44 Ca/ 40 Ca 
positive shifts in some samples going up to +181 (see their Figure 5). They fully recognized 
the possible role of AGB processing and showed the pattern to be expected for a 2 M 0 star. 
Our calculations of both pattern and yield are in accord with their report. They did not 
include 46 Ca and 48 Ca, but they recognized that 46 Ca would be produced. They further 
discussed the relationship between 42 Ca/ 40 Ca and 43 Ca/ 40 Ca to be expected from AGB 
sources. 


In considering the problem of 48 Ca we also recognize that any inventory of this isotope 
must be due to some SNe sources. However, the existence of both positive and negative 
effects for 48 Ca (and many other isotopes) proves that any simple addition mechanism 
involving one isotope is not possible. Inspection of the curves of Q it40 (Figure [5]) or of E lW 
(Figure [6]) shows that no 48 Ca is produced in AGB processing and that this nucleus is also 
not efficiently destroyed. We further note that the production of 42 Ca, 43 Ca, and 44 Ca is 
subdued, but is a natural result of AGB processing. For the case of a Solar initial state 
in a star, the results enhance 42 Ca relative to 44 Ca. This can be clearly seen in Figure [9j 
where it is shown that enrichment for 42 Ca is always greater than for 44 Ca. It follows that 
shifts in 42 Ca/ 44 Ca are to be expected. As all of the estimates of 48 Ca excesses/deficits from 
measurements use a normalization to Solar 42 Ca/ 44 Ca, it is clear that any enhancement of 
42 Ca/ 44 Ca by neutron captures will result in a calculated excess of 48 Ca. Further, if there 
are small variations 42 Ca/ 44 Ca in the ISM at the formation of a star, this will in general 
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result in either excesses or deficits in the calculated 48 Ca/ 44 Ca ratio 


From the observational data, we note that all Ca isotopic effects are small, being 
typically a few e units in bulk meteorites and most CAIs. In oxide CIRCONs, the effects 
are less than 100 per mil in most samples. As discussed above, sev eral of the cases show a 


rather strong resemblance to the AGB model (see also 


Nittler 


2005) and may be explained 


by variability in the initial inventories of the AGB stellar sources as argued for Ti. 


In considering the Ca problem we focus on three issues: 

1) There is no production and no effective destruction of 48 Ca in AGB stars; 

2) There is a high degree of enrichment in 46 Ca as a result of AGB processing and very 
few measurements of this rare isotope. The only clear case is the data on EK141 with only 
modest excesses of 46 Ca (3 per mil) with a much larger excess in 18 Ca (13 per mil); 

3) The effect of s-processing on Ca is to produce modest enrichments (see Figures [5] and 
m i 11 42 Ca, 43 Ca, and 44 Ca with the net effect in all cases to have larger enrichments in 42 Ca 
relative to 44 Ca. 


If the Ca calculation is correct, then there should be very large enrichments in 46 Ca. 
The large effects we predict for 46 Ca are also to be found in the FRUITY data base. This 
isotope has only rarely been measured, but no positive evidence of large enrichments in 
46 Ca are even hinted at. If the lifetime of 45 Ca is not the laboratory value, with t\/ 2 = 163 
days, but is controlled by short li ved isomeric states, or by other processes affecting weak 


interactions in stars flSimonucci et al. 


2Q13|), then the calculated large effects on 46 Ca and 


the inferred large shifts are not valid. If this is not the case, then we are obliged to infer 
that there are very large 46 Ca effects to be found or that the s-processing fails to explain 
the Ca data, which would place in jeopardy the approach presented here of processing of 
Fe group elements in AGB stars as the cause of observed shifts in these elements. The 
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behaviour of 48 Ca is a matter of some issue. The measurements of 48 Ca in almost all cases 
depend on a normalization procedure using 42 Ca and 44 Ca. In so far as the material being 
analyzed has small variability in the Ca isotopes, then the normalization 42 Ca/ 44 Ca must 
remain in doubt. If a sample has been subject to processing in AGB stars, then it will be 
enhanced in 42 Ca relative to 44 Ca. A 2 e units enhancement in the 42 Ca/ 44 Ca ratio in such 
material would give a calculated 4 e units excess in 48 Ca/ 44 Ca. We consider that this may 
be an explanation of some of the reported 48 Ca excesses. However, as pointed out by D.A. 
Papanastassiou (private communication) if the shifts in 42 Ca/ 44 Ca are due to s-processing, 
then there will be distinct effects on 43 Ca that are not found. Some variability in the 
48 Ca inventory in the ISM is also a most plausible scenario; however, the issue of the 48 Ca 
anomaly being due to a normalization calculation procedure requires serious attention. 

We find the disparity between the patterns observed for Ca and the one to be expected 
from AGB processing of an assumed initial solar composition to be similar to what was 
cited above for Ti. As the shifts in 12 Ca , 43 Ca, and 44 Ca are small, due to the low neutron 
capture cross section of 40 Ca, small variations in the initial stellar composition will have 
much enhanced effects on the patterns observed. The Ca problem was one of the most 
serious difficulties we had in addressing the role of AGB re-processing. 


4.7. Zn 


The only data on Zn are for CAIs, reported by 


Voelken ing & Pa panastassiou (119901 ). 


Normalizing the data to 64 Zn/ 68 Zn, they found an excess in 66 Zn of 16 e units in EK141 and 
deficits of ~ 10 — 20e units for 70 Zn, with large errors. 67 Zn was not measured for technical 
reasons. All other CAIs showed normal Zn. They discussed the extent to which these 
anomalies could be produced in neutron rich NSE synthesis. As 68 Zn is the highest enriched 
Zn isotope from s-processing (see Figures [5] and [6]) and all the Zn isotopes are significantly 





enriched by s-processing, any normalization procedure is complex. In the absence of data 
for 67 Zn in EK-141, a relationship to the s-processing model is made difficult. No clear 
explanation for these results are yet in hand. 


5. Summary and conclusions 

This work was initially directed to establish whether the excesses observed in 54 Cr, 
58 Fe, and 64 Ni in some CAIs could be the result of AGB processing of matter from the 
ISM. Detailed calculations showed that the data could be the result of neutron captures in 
such stars and further, that the effects on these isotopes were much larger than on any of 
the other isotopes of these elements. As a result, the large excesses in 54 Cr, 58 Fe, and 64 Ni 
would be present even if there were substantial variations in the initial isotopic composition 
of the contributing AGB stars. These “Fe group” elements were originally made in SNe, but 
their products have been subsequently reprocessed in AGB stars, resulting in the observed 
effects. Assuming this to be generally applicable to other “Fe-group” elements, we then 
sought to establish whether this might explain the observations on Ti, Si, and Ca. We 
therefore expanded the presentation to include the wealth of data on pre-solar circumstcllar 
grains. 

For Ti, the situation is roughly similar to that for Cr, Fe, and Ni, where the 
predominant enrichment is in one isotope. However, in addition to large effects predicted 
for 50 Ti, there is also significant to major production of 46,47,49 Ti relative to 48 Ti. The large 
effects of 50 Ti produced were used as a guide. Using the important and extensive results 
on circumstcllar condensates (CIRCONs), known from other workers to be associated with 
AGB sources, we sought to relate the Ti isotopic data to AGB processing. It was shown 
that a substantial part of the Ti data on CIRCONs show clear effects of s-processing in stars 
having close to Solar initial Ti isotopic ratios. However, as is well known, many variations 
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in the Ti data did not fit into the AGB enrichment pattern, particularly for grains that 
were not greatly enhanced in 50 Ti. Previous workers had inferred that the scatter in Ti 
isotopic abundance patterns in CIRCONs must be due to variations in the initial states of 
the different stellar sources. We are in accord with this, but are of the view that, in general, 
the Ti data are the result of very substantial AGB processing, superimposed on a base of 
a variable initial isotopic inventory from the local ISM from which the star formed. The 
required variability is in all cases small: ~ 10% in isotopes that are at most ~ 10% of the 
total Ti inventory (for example, 46 Ti changing from an 8% to an 8.8% of total Ti in the 
local ISM). The negative effects observed can be explained by an increase in 18 Ti relative 
to the other isotopes by less than 1%. 


With regard to the widely recognized 48 Ca anomalies, commonly attributed to SNe, 
we consider that the assumption of Solar 42 Ca/ 44 Ca ratio in normalizing the data could 
generate apparent 48 Ca / 44 Ca excesses/deficiencies. An AGB model predicts large excesses 
of 46 Ca which are not observed. A possible explanation may result from the effects of the 


20131), or from a decay 


high stellar electron densities on weak interactions (jSi monucci et ah 
rate of 45 Ca governed by some short lived isomeric state in stellar conditions. Work to 
clarify these possibilities is in progress. 


For Si from AGB sources, the regular trends of £ 30,28 versus $ 29,28 with a slope of ~ 1 
found for the “main stream data”, do not appear to reflect the effects of neutron captures 
with the standard cross sections or neutron exposures. As noted above, the Si isotopic 
shifts to be expected from AGB processing are very small compared to those for Ti. Any 
heterogeneity in the ISM such as found for Ti would overwhelm any Ca shifts, so these 
two data sets should be considered decoupled. The observed Si-Ti isotopic correlation has 
earlier been considered by many authors to be the result of Galactic Chemical Evolution. 
However, considering the heterogeneity in initial states clearly shown by the Ti results and 
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the wide range of Si isotopic abundances from different SNe sources (particularly SNela) 
and the small Si effects from AGBs, we infer that the slope ~ 1 is simply due to a somewhat 
variable (~ 10 to 20%) contribution of SNela to the various local regions where these stars 
formed. Considering the fact that from 1/2 to 2/3 of all the Fe nuclei must come from these 
sources, which produce 28 Si and very little 29 Si and 30 Si, this a ppears to be a plausible 


cause of the observed trend, as early suggested by 


Lugaro et ah 


(119991 ). We consider the 


granularity of the ISM composition to be the greater governing agent for most of the data 
reported on Si and Ca and any coupling between s-processing for Ti, Cr, Fe, and Ni with 
Si and Ca will be weak, because of the low effects resulting from s-processing on silicon and 
calcium themselves. 


We note that the calculations, both in the stellar models and in the data representations, 
use Solar reference isotopic values. As derived from equations (9), (12), and (13) the 
enrichment factor relative to the solar reference (expressed by the factor ay*,) is (see 
equation (16): 


(10- 3 6X ~ <*«) / (10 - 3 Sf,t - a jit ) = (Q® - o«) / (Qf :t - a M ) 


This is calculated for a sample formed in the envelope of a star that did not have an initially 
solar isotopic composition (see section 2.3.3 for details). This simple relation connects 
the calculated laboratory data relative to Solar standards to the value that would be in 
a model star of solar initial isotopic admixture, but is actually offset by the actual initial 
composition through the abundance factor ay*,. This permits of determining the initial 
composition of the source if one ay*, can be estimated. A shift due to ay k will be very 
significant (e.g. for ay fc = 10 -2 , this is a 1%, or 10 2 e unit shift in the isotopic ratio). 


Support for possible isotopic variability may be fo und in the known variability in 


chemical abundances for stars with [Fe/H] ~ 0 (see e.g. 


Majorca et ah 


2012 


Martig et ah 


20151 ). Certainly, 10% to 20% shifts must be expected in the ISM. The case for galactic 











39 


changes will be more clearly defined when the observational data and theory require 
metallicities much below solar and hence significantly earlier epo chs. This may be 


applicable for some of the oxygen results in the report by 


Niftier & Coswik (1997). The 


role of heterogeneity that we invoke has been widely known, but only in connection with 


much larger scale issues of Galactic Chemical Evolution. As noted by 


Niftier et ah 


(2QQ8) 


(section 4.2.3): “The fact that the total spread in the grain data is larger than is expected 
for nucleosynthesis in a single O-rich AGB star suggests that much of the range indeed 
reflects variations in initial compositions du e to GC E. This is anal 


and Ti isotopes in presolar SiC grains (1Alexander & Niftier 


is ana 

Jogous to the case 

1999; 

Zinner et al. 

2007) 


concur that the variations reveal variable initial compositions; however we believe that they 
reflect the intrinsic granularity and heterogeneity of the local ISM from which stars form 
and, to a lesser extent, some effects of broader-scale Galactic Chemical Evolution. For this 
last, the issue is to identify CIRCON grains which require low metallicity to explain their 
isotopic data. 


The ISM does not have an exactly solar isotopic composition, but is heterogeneous 


both in time and space. That there is a difference in the initia' 


has long been recognized by many workers flAmari et ah 


2001 


composition of the sources 


2004; 


Huss fe Smith 


2007). 


We agree with previous researchers that the silicon and titanium isotope arrays exhibited by 
presolar SiC grains reflect primarily the initial co mpositions of the parent stars, overprint ed 


by a component of s-process nucleosynthesis (see 


Arnari et al. 


2001 


Huss fe Smith 


2007L 


However the small s-components are the effect of early stages of dredge up, which emphasize 
the inherited component; while more advanced stages are dominated by the s-process itself. 


Measurements of the isotopes 54 Cr, 58 Fe, and 64 Ni are very difficult using SIMs 
techniques, because of the isobaric interferences. However, some aspects of the AGB effects 
can be reconstructed with complete data sets of Ti, Cr, Fe, and Ni on single samples of 
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CAIs. Insofar as the isotopic effects are sufficiently large, it should be possible to establish 
whether the predicted isotopic correlations are really found. For example, the absence of 
58 Fe in the presence of excess 64 Ni or of high 50 Ti values would cast serious doubts on the 
approach presented here. 

A careful assessment should now be made of the codes for s-process in AGB stars with 
a focus on the Fe group elements, so that errors in the nuclear physics or the codes are not 
propagated in our models. More broadly, the general consequences for all nuclides of the 
results of re-processing ISM matter in AGB stars requires careful attention for all affected 
nuclei, including radioactive isotopes. Attempts at identifying contributions from Galactic 
Chemical Evolution (not from local granularity) should be focused on CIRCONs grains 
that may reasonably be considered to come from stars with [Fe/H] values well below zero. 
Very high enrichments of 50 Ti may aid in identifying such grains. Further work in this held 
may also be directed to look for correlations between some of the signatures discussed in 
the present paper and the several recent reports on very small isotopic variations in many 
heavier elements (eg Mo, Ru, Hf, and W), which show s-process excesses or deficits in 
leaches of various meteoritic materials. The problems associated with leaching experiments 
for different chemical elements will of course remain. 
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Figure Captions 


Figure [0 The network of reactions we considered to be active for the iron group 
nuclei relevant for our discussion. In the plot we show approximate 30keV values for 
cross sections (lower rows in the boxes indicating the nuclei) expressed in millibarns and 
laboratory lifetimes for unstable nuclei (upper rows). Note that these data are indicative: 
cross sections vary in a complex way as a function of temperature and also some lifetimes of 
radioactive species do the same. Note also that our network is, to our knowledge, exactly 
the same as adopted by the FRUITY database quoted in the text. Note also that our 
network is, to our knowledge, the same as adopted by the FRUITY database . 


Figure [2[ The enrichment factors (Ef k ) for Cr, Fe, and Ni for the envelopes of model 
AGB stars (with mass of 1.5 and 3M 0 ) having initially both a^Solar metallicity, Z = Z 0 , 


and Solar isotopic ratios, following the models by 


Trippella et al. 


(2014), but with the 


network updated for the Fe group elements as given in this study. Note that the calculation 
of 54 Cr and 58 Fe from our report lies well below that found in FRUITY for 1.5M 0 . Figure 
[2] shows the cases for 1.5 and 3 M r . : for different pulses and includes some results from the 

(201l|, using their estimates for the enrichment 


FRUITY on-line database 


Cristallo et al. 


factors in the envelope. The right part of the figure (Figure El Panels d, e, and f) shows, 
for clear comparison, the enrichment factors for each case normalized to E® k =1 where 
“j” represents the isotope with the largest enhancement ( 54 Cr, 58 Fe, and 64 Ni). In our 
calculation the major shifts in isotopic abundance are clear. There is good agreement in the 
isotopic patterns for Cr and Fe for the results presented here and those from the FRUITY 
data base. There is a discrepancy in the case of Ni, due to the underproduction of 64 Ni 
obtained from the FRUITY data base. 
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Figure [3] The enrichment factors (Ef k ) for the envelopes of model AGB stars with 
initial Solar isotopic ratios, but for Z = Z Q /3. See details in the caption of Figure [21 Note 
that here, the data points for 58 Fe and 54 Cr in the 1.5M 0 case by both calculations are 
close. 


Figure [4j The F,^.-values for Ti at the last pulse from the FRUITY outputs, for 
two values of the metallicity that bracket our choice (Z 0 = 0.0154). As a comparison, the 


results we o 
Figure 1 in 


atained in this p aper , both for our standard choice of the 13 C pocket mass (see 


Trippella et al. 


2 0141) and for the smaller pocket of 0.001 M & . The isotopic 


patterns for both calculations are essentially the same. However the much larger effects for 
the Fruity output as compared to our calculations must be investigated. 


Figure [5] Isotopic enrichments for a wide suite of elements of s-processed matter in a 
parcel before mixing in to the envelope Qf k (see section 2.3.1 in text). The values of Qf k 
are arbitrarily chosen here at the last pulse in a 1.5 M 0 model. Solar initial isotopic values 
were used in all these calculations. The index isotopes u k” used are 40 Ca, 48 Ti, 52 Cr, 56 Fe, 
58 Ni, and 64 Zn. Comparison should be made between the patters represented by Qj^ and 
the envelope enrichment patterns shown in figure [6l It is important to recognise that to 
reach around unity would require, from our calculations, a more massive star and low 

metallicity (Z 0 /3) source (see numerical numbers in Figure [6]). Full red triangles represent 


the Q values for smaller 1 3 C pocket 


(2014). 


in or der to reproduce the effects suggested by e.g. 


Gallino et al. 

1998; 

Wasserbure et al. 

2006 

) from the Case A calculations of 

TrioDclla et al. 
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Figure [6j Summary of the isotopic compositions for a wide variety of “Fe 
group” elements in the envelope of AGB stars. We have represented the results normalized 
to F® fc =1 for those isotopes “j” that are most enhanced by s-processing in our calculations. 
These are, respectively, 46 Ca, 50 Ti, 54 Cr, 58 Fe, 64 Ni, and 68 Zn. The index isotopes “fc” used 
are 40 Ca, 48 Ti, 52 Cr, 56 Fe, 58 Ni, and 64 Zn. Notice that the form of these curves is almost 
independent of the pulse, the stellar mass or the metallicity. The curves shown here should 
be compared with the Qf k shown in Figure O Ca shows the most discrepant behaviour. 
Otherwise, the curves are almost indistinguishable. See text for discussion of the individual 
elements. Note the peak at 46 Ca, resulting from using the laboratory lifetime of 45 Ca. This 
effect is not supported by observations. 


Figure [7] Measurements on Ti with quoted errors (1 Amari et ah 


2Q01) for the CIRCON 


grain KJGM1-158-5, which has the largest measured enrichment in 50 Ti. For comparison, 
we show the curve for Qf 48 in star a of 1.5M 0 and Z Q , which is normalized to £ 50 ,48 = 
990 per mil. Note the very good agreement between these data and AGB s-processing. 
It’s important to recognize that to reach F 50 ,48 around unity would require, from our 
calculations, a more massive (3M 0 ) star and low metallicity (Z Q / 3) source. 


Figure [8] A wide selection of Ti isotopic data on main stream CIRCONS compared 


to A GB s process patterns. These samples 


from 


Amari et al. 


tail), 


Amari et al. 


lave 5k down to 200 per mil. Data taken 


(120041) and 


Huss & Smith (12007! ). For comparison 


we show the Qi^s normalized to Q 50,48 = 1 (that is Qi,A&/Qho, 4 s)- The data shown are in the 
form (£ i)48 ± e i!48 ) / (£ 50 , 48 ) where e^g is the uncertainty quoted for that ratio. No errors 
were propagated from the values for £ 50 , 48 - It can be seen that the data cluster around 
the AGB model calculations, but with some clear differences. In several cases they are 
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well within the error envelope, but there are also many exceptions. For example, sample 
number 320 has clear excesses of 46 Ti and 4 ' Ti relative to the model, but matches 49 Ti very 
well. Sample KJGM4C-314-3 appears to be quite discrepant, with a large excess of 49 Ti. 
However, if one considers the errors in both 50 Ti/ 48 Ti and 49 Ti/ 48 Ti, the results could be 
compatible with AGB processing. 


Figure [9l The enrichment factors for the Ca isotopes 42 Ca, 43 Ca, and 44 Ca, assuming a 
solar initial isotopic composition. 46 Ca and 48 Ca are omitted. These results are normalized 
to Ef 240 = 1; otherwise the display is same as in Figure [HI It can be seen that Ef 240 and 
Ef 3 40 are enhanced relative to Ef 4 40 in every case. Assuming 42 Ca/ 44 Ca as normal, then a 
calculation of the 48 Ca/ 44 Ca ratio will result in an apparent excess of 48 Ca for s-processed 
material; any other non-Solar source will give effects in 48 Ca/ 40 Ca. It is seen that these 
results for Ca are well behaved (compare to Ca in Figure [6]). 
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Fig. 4. 
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Fig. 7. 
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